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Abstract 

Starting from tiie fluctuating Boltzmann equation for smooth inelastic hard spheres or disks, 
closed equations for the fluctuating hydrodynamic fields to Navier-Stokes order are derived. This 
requires to derive constitutive relations for both the fluctuating fluxes and the correlations of the 
random forces. The former are identified as having the same form as the macroscopic average 
fluxes and involving the same transport coefficients. On the other hand, the random force terms 
exhibit two peculiarities as compared with their elastic limit for molecular systems. Firstly, they 
are not white, but have some finite relaxation time. Secondly, their amplitude is not determined 
by the macroscopic transport coefficients, but involves new coefficients. 

PACS numbers: 45.70.-n,05.20.Dd, 05.60.-k,51.10.+y 
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I. INTRODUCTION 



Granular matter is ubiquitous in Nature. It has a tremendous impact on a wide range 
of industries and also raises important conceptual challenges. When the grains move freely 
and independently between collisions the system is referred to as a granular gas The 
simplest model of a granular gas at the particle level of description is an ensemble of in- 
elastic hard spheres or disks |2|. In the last decade or so, a great deal of progress and 
understanding of granular gases has been achieved ^, 4| by using this model, specially in 
the dilute limit in which the average behavior of the gas is accurately described by the inelas- 
tic Boltzmann equation js, 6|. This includes derivation of the hydrodynamic Navier-Stokes 
equations with explicit expressions for the transport coefficients, study of several Newto- 
nian and non-Newtonian steady states, analysis of the spectrum of the linearized inelastic 
Boltzmann operator, identification of several instabilities both in freely evolving and heated 
granular gases, and detailed characterization of the distribution function of several relevant 
states The theoretical predictions have been corroborated by molecular dynamics and 
direct Monte Carlo simulation results. On the other hand, knowledge about fiuctuations 
and correlations in dilute granular gases is much more limited, although they are known to 
play a fundamental role in the macroscopic behavior of granular fiows in many cases. 

The purpose of this paper is to derive Langevin-like equations describing the dynamics 
of the fiuctuating hydrodynamic fields in a freely evolving homogenous dilute granular gas, 
namely in the usually termed homogeneous cooling state (HCS) [Sj. Of course, in the elastic 

jydrodynamic equations proposed by Landau 



|. The starting point will be the fiuctuating 
I by using non-equilibrium statistical me- 



limit the equations reduce to the fiuctuating 
and Lifshitz for equilibrium molecular fluids 
Langevin-Boltzmann equation derived in ref. 
chanics methods. The equation has the same mathematical form as the linearized inelastic 
Boltzmann equation plus and additive white noise term. 

Hydrodynamic fiuctuations in freely evolving granular gases have been already partially 
investigated. First studies used mesoscopic equations in which the Landau and Lifschitz 
fiuctuating hydrodynamic equations for molecular gases were directly modified to incorpo- 
rate the effect of energy dissipation in colhsions, the resulting equations being expected to 



be valid in the quasi-elastic limit 



In particular, the correlations of the noise terms were 



assumed to be proportional to delta time functions (i.e. white noises) and determined by the 



2 



Navier-Stokes transport coefficients, using the same expressions as for molecular systems. 
Also, some results based on a single relaxation model kinetic equation have been reported 



12| . In the above studies the interest focused on the initial buildup of spatial correlations of 



the hydrodynamic fields near the so-called clustering instability. Approaches starting from a 
more fundamental description of the granular gas, i.e. the Liouville equation of the system, 
have also been carried out. In particular, fluctuations and correlations of the total energy 
of an isolated granular gas have been computed and a good agreement between theory and 
simulations has been found [3, Q • Moreover, in ref . jioj] , the fluctuations and correlations 
of the transversal component of the velocity field were studied in detail. It was shown that 
the Langevin equation obeyed by this field in a granular gas has two crucial differences 
with the elastic limit. Firstly, the noise is not white but has a finite correlation time and, 
secondly, the amplitude of its second moment is not determined by the shear viscosity but 
by some new coefficient. It must be noticed that the transversal velocity field is special in 
the sense that it evolves uncoupled from the other hydrodynamic fields. 

The purpose of this paper is to extend the analysis in p^], by deriving the whole closed 
set of fluctuating hydrodynamic equations for a dilute granular gas in the HCS. In this 
context, it must be realized that this state is the homogeneous reference state for granular 
gases, playing a role somehow similar to the equilibrium state in molecular fluids. Therefore, 
understanding fluctuations and correlations in this state is an unavoidable first step towards 
a deep knowledge of them in more complex and also more realistic states of granular gases. 

It is worth to mention also some previous studies of fluctuations dealing with driven 
granular gas models, in which the grains are assumed to be in contact with some external 



energy source or thermal bath 15Nl7l| . Since the latter modifies the stochastic properties of 



the granular gas in a nontrivial way, it is not evident that there is a direct relation between 
fluctuations in driven models and the free model being considered here. 

The plan of the remaining of the paper is as follows. In Sec. HIl the Boltzmann-Langevin 
equation for smooth inelastic hard spheres or disks is shortly reviewed, and the appropri- 
ate dimensionless time and length scales are introduced. From the above kinetic equation, 
balance equations for the fluctuating hydrodynamic fields in the HCS are directly derived, 
as shown in Sec. IIIII These equations are worthless until they are closed by deriving con- 
stitutive relations for the fluxes appearing in them and also explicit expressions for the 
correlations of the noise terms. This is done here by introducing a projection operator over 
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the hydrodynamic subspace of the distribution functions. This subspace is generated by the 
hydrodynamic eigenf unctions of the hnear inelastic Boltzmann operator, whose definitions 
and expressions are also briefly reminded. The equation for the number density involves 
neither flux nor noise term. The Langevin equation for the velocity field is derived in Sec. 
IIVI and the one for the energy field in Sec. El In both cases, expressions for the fluctuat- 
ing fluxes and for the correlation functions of the fluctuating forces are provided. These 
expressions are not exact, but have been obtained under well defined and controlled approx- 
imations. It is shown that the fluctuating forces have finite relaxation times which can be 
related with non-hydrodynamic eigenvalues of the linearized Boltzmann operator. For the 
sake of clarity, many details of the calculations as well as technical points are given in six 
appendixes. Finally, Sec. I VI I contains a short summary and some final remarks. 



II. BOLTZMANN-LANGEVIN EQUATION FOR SMOOTH INELASTIC HARD 
SPHERES OR DISKS 

The system being considered is a dilute gas of smooth inelastic hard spheres {d = 3) or 
disks {d = 2) of mass m and diameter a. Inelasticity of collisions is modeled by means of 
a constant, velocity independent, coefficient of normal restitution a, defined in the interval 
< a < 1. At a mesoscopic level, the system is described by a fluctuating or stochastic one- 
particle distribution function, F{r,v,t), whose average, denoted by f{r,v,t), is the usual 
one-particle distribution function giving the average number of particles density being at 
position r with velocity v at time t. This average distribution obeys the nonlinear inelastic 
Boltzmann equation 5|. 

The system is assumed to be in the homogeneous cooling state (HCS), macroscopically 
characterized by a uniform number of particles density n, a vanishing flow velocity, and a 
uniform granular temperature Tuit) that decreases monotonically in time due to the energy 
dissipation in collisions. To describe fluctuations about the HCS, it is convenient to introduce 
dimensionless length and time scales defined by 

V 

and 

(2) 
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respectively. In the above expressions, A = {na'^ ^) ^ and 



vo{t) 



m 



(3) 



Then A is proportional to the mean free path of the gas in the HCS and vo{t) is a characteristic 
thermal velocity. The granular temperature is defined from the average kinetic energy of 
the grains with the Boltzmann constant formally set equal to unity. It is easily verified that 
the time scale s is proportional to the accumulated average number of collisions per particle 
occurring in the system in the time interval between and t. Consistently with Eqs. ([1]) 
and (E]) a new velocity scale is introduced through 



V 



voit) 



Also, the deviation of F{r,v,t) from its average in the HCS, fnivjt), 

6F{r,v,t)=F{r,v,t)-fH{v,t), 



is transformed to a dimensionless form by defining 



6F{£, c, s) = n-^v^{t)6F{r, v, t) 



This function obeys the Boltzmann-Langevin equation lOj 

6F{£, c, s) = S{£, c, s) 



9 5 w X 



(4) 



(5) 



(6) 



(7) 



where A(c) is the linearized inelastic Boltzmann operator given in Appendix |A] and the term 
S{£, c, s) has the properties of a white noise term. 



for s > s', and 



{Si£,c,s))H = 0, 
{S{£,c, s)6F{£',c',s'))h = 0, 



{S{£, c, s)S{£', c', s'))h = n-^\-'^5{£ - £')6{s - s')r{c, c' 



(8) 
(9) 

(10) 



Here and in the following, angular brackets with the subindex H are used to denote stochastic 
average in the HCS. The amplitude r(c, c') of the noise is given by 



r(c, c') = - [A(c) + A(c')] x{c)5{c - d) + To(c, c')x(c)x(c'), 



(11) 



with Tq[c,c') being the inelastic binary colhsion operator given in Eq. ( lAip and x(c) the 
dimensionless scaled velocity distribution of the HCS defined as 

x{c)^n-X{t)fH{v,t). (12) 

An accurate expression for this distribution is known and it is reminded in Appendix [Xj 
Eq. flA6|) . For the present purposes, it is useful to employ the Fourier representation and to 
introduce 

SF{k,c,s) = J de e-'^-^5F{l, c, s) . (13) 
In Fourier space, Eq. ([7]) becomes 

5F(fc,c,s) = ^(fc,c,s), (14) 
where 

A(fc,c) = A(c) -ifc • c (15) 

is the linear inhomogeneous inelastic Boltzmann operator. The property of the noise in Eq. 
( ITOjl now reads 

{S{k, c, s)S{k\ c', s'))h = ^5k,-k'5{s - s')f (c, c'), (16) 

^k-k' being the Kronecker delta symbol and V = VX^''' the volume of the system in the 
length scale defined by A. 



gj-A(fc.c) 



III. BALANCE EQUATIONS FOR THE FLUCTUATING HYDRODYNAMIC 
FIELDS 

Consider the mesoscopic number of particles density N{r, t), momentum density G{r, t), 
and energy density E{r,t). Dimensionless deviations from their average values in the HCS 
are introduced through the definitions 

Sp(k, s) = ^^^'""^ = [ dcSF{k, c, s), (17) 
n J 

Suj{k,s) = ^^ihlA = f dccSF{k,c,s), (18) 
mnvo{t) J 

Hk, s) = ^^1^ [E{k, s) - ^riTnis)] = -^j dcc^6F{k, c, s). (19) 
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In [10|, balance equations for these fields were derived by taking velocity moments in the 
fluctuating Boltzmann equation. In Fourier space they read 

d 

— Sp{k,s) + ik-5uj{k,s) = 0, (20) 

- I) '^'^(^' ^) + • s) = 0, (21) 

^ 6e(k, s) + i k ■ 6u;(k, s) + 5C(k, s) + i-k- S(f)(k, s) = SJk, s). (22) 
OS a a 

In the above equations, Co is the cooling rate of the HCS. Its definition is given in Appendix 
lAt where an approximated expression is also provided. Moreover, 6r\{k,s) and (5</)(fc, s) 
are the fluctuating pressure tensor and heat flux, respectively. They are functionals of the 
fluctuating distribution function, 

6n{k, s) = + j dcA{c)6F{k, c, s), (23) 

S(f){k,s) = J dk'E{c)5F{k,c,s), (24) 
where I is the unit tensor of dimension d and 

A(c) = cc-^l, (25) 



2 d + 2 



S(c) ^i^c^-^jc. (26) 

The equation for the fluctuating energy field, Eq. f l22|) . involves two terms vanishing in the 
elastic limit a — )■ 1. One is associate with the cooling rate fluctuations 

s) = y dcc^A{c)SF{k, c, s), (27) 

and the other one is the noise term Se{k, s) arising directly from the fluctuations in phase 
space, 

S,{k,s) = dcc^S{k,c,s). (28) 
From Eqs. ([8]) and ( 1T6|) it follows that 

{S,{k,s))n = (29) 

and _ 

{S,{k, s)S,{k', s'))h = ■^Sk,-k'S{s - s') I del dc'c\'^r{c, c'). (30) 



Also it is 



and 



{S,{k,s)F{k',c',s'))H = (31) 



(5,(fc, s)Sp{k', s'))h = {S^ik, s)5Lj{k', s'))h = {S^ik, s)5e(fc', s'))h = 0, (32) 



for s > s'. 

Equations (!20|) - (!22|) are not closed, since they contain the quantities 511, 6<j), and 6(, 
defined above in terms of 6F{k, c, s), as well as the noise term Se{k,s). The aim in the 
following will be to obtain a self-consistent description for the fiuctuating fields. In this 
context, it will be useful to consider the eigenvalue problem for the inelastic homogeneous 
linear Boltzmann operator 18|, ll9|] , 



(33) 



The solutions of this equation corresponding to the infinite wave length limit (A; = 0) of the 



hydrodynamic equations are given by fisl. Il9| 



Ai = 0, A2 = | 



A.-? — — , 



(34) 



d 



ei(c) = x(c) + ^ ■ [cx(c)] , 12(c) 



e3(c) = -^-Mc)]. (35) 



The eigenvalue A2 is rf-fold degenerate. Normalization of the distribution function re- 
quires that the velocity distribution functions being considered be integrable. However, 
the strongest requirement is made now that they be elements of a Hilbert space with scalar 
product defined as 

{g\h) ^ [ dcx-\c)g*{c)h{c), (36) 



where g*{c) is the complex conjugate of g{c). The operator A(c) is not Hermitian and the 
eigenf unctions ^/3(c) are not orthogonal. Then, it is convenient to introduce a set of functions 
,^^(c) biorthogonal to the above eigenfunctions, i.e. verifying 



(37) 



A convenient choice is 



18 



19| 



ei(c)=x(c), Uc) = cx{c), e3(c) = ( ^ + ^ ) X(c) 



(38) 
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Using the bi-orthogonal sets of functions, a projection operator V over the hydrodynamic 
part of the Hilbert space can be defined by 

<i+2 

P^(c)^ J]e^(c)(e^|^7). (39) 

/3=1 

By means of P, the fluctuating one-particle distribution function can be decomposed into 
its hydrodynamic and non-hydrodynamic components, 

5F{k, c, s) = V6F{k, c, s) + V±6F{k, c, s), (40) 

where V± = 1 — V. Application of V± to both sides of Eq. ( IT^ and formal integration of 
the resulting equation yields 

V±SF{k, c, s) = U{k, c, s)V±6F{k, c, 0) + [ ds'U{k, c, s')V± \-ik ■ cVSF{k, c, s - s') 

Jo 

+S{k, c, s - s')j , (41) 

where 

W(fc, c, s) = exp [sV±Aik, c)Vi_] . (42) 

The hypothesis is made now that for large enough s, the first term on the right hand side 
of Eq. f l4T|) becomes negligible. This is related with the ageing to hydrodynamics, implying 
that the non-hydrodynamic part of the initial condition is forgotten on the hydrodynamic 
time scale. In addition, the limit of small wavevector k is considered and only terms up to 
first order in it are kept. Then, for large values of s, Eq. (141 p becomes 

V±6F{k,c,s) - [ ds'T^e''^^""^ {-ik ■ c)V6F{k,c,s - s') 
Jo 

+ [ ds'U{k,c,s')V±S{k,c,s- s'), 
Jo 

(43) 

This expression will be used in the next sections to obtain explicit expression for the fluc- 
tuating fluxes and for the cooling rate in the Navier-Stokes approximation. 

IV. LANGEVIN EQUATION FOR THE VELOCITY FIELD 

By direct calculation, it is easily verified that 

J dcA{c)C0{c) = O. (44) 
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Here ^/3(c) stands for any of the hydro dynamic eigenfunctions of A(c) given in Eq. (135|) . 
Therefore, 

^ dcA{c)VSF{k,c,s) = (45) 



and, consequently, 5F can be substituted by V±6F on the right hand side of the expression 
for the fluctuating pressure tensor, Eq. (123!) . Using next Eq. (HSj) gives 



where 



and 



su{k, s) ~ I + 5^n(k, s) + R(k, s), 



5in(fc, s) = ds j dcA{c)e''^^''\-ik ■ c)V5F{k, c,s- s') 



R(k,s)= I ds' I dcA{c)l({k,c,s')V±S{k,c,s-s'). 



Equation fl2T|) implies that 



Suj{k,s-s') ~ e-'''^°/^Suj{k,s) 



(46) 

(47) 
(48) 

(49) 



to lowest order in k. Using this, it is obtained that 10| 



6illij{k,s) = -ir]{s) 



ki6ujj{k, s) + kj6ui{k, ~ ^ij^ " S^{k, s) 



(50) 



valid to first order in k. In the above expression, ri{s) is the time-dependent dimensionless 
shear viscosity of the HCS |20|, 

d 



') = + j dcAiJ{c)<^2,^Jic, s), 

$2,i(c,s)= f ds' e^'^^^-^^'^^^Af^yy 

Jo 



(51) 



(52) 



Equation fl5T|) agrees with the result found from the nonlinear Boltzmann equation by using 
the Chapman- Enskog algorithm, and its long time limit has been evaluated in the first 



Sonine approximation 



21 



221. 



The last term on the right hand side of Eq. (1461) is defined in Eq. (1481) and it has the 

property 

{R{k,s))H = 0, (53) 
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that follows from Eq. (|8]). Moreover, in Appendix [B] it is shown that for s S> 1, s' ^ 1 and 
to lowest order in /c, it is 

, V"^ , 2 

{Rij{k, s)Rim{k', s'))h = —Sk-k'G{\s - s'\){6uSj,m + SimSjl - -SijSim), (54) 

where 

G{s) = ^^-L_J2 [ dc^ f dc2 A,,(ci)A,,(c2)e^^(^^)0^^(ci,C2), (55) 

with 0j^(ci, C2) being the solution of the equation 

[A(ci) + A(C2)] M^u C2) = -Pf f (ci, C2). (56) 

The operators V^^^ and V^^^ are defined like V±, but acting on functions of the velocities Ci 
and C2, respectively. An approximation will be introduced at this point. The noise amplitude 
r(ci, C2) given by Eq. (ITT]) has two contributions of a rather different physical origin. The 
first one reflects fluctuations induced by collisions as a consequence of two different particles 
colliding independently at the same position with the same environment. This contribution 
does not vanishes even for a molecular gas at equilibrium 2^ . On the other hand, the second 
contribution to r(ci,C2) is directly related with the velocity correlations between particles 
in the HCS and vanishes for an ordinary fluid at equilibrium. Here, the hypothesis will be 
made that only the hydrodynamic component of the velocity correlations, as extracted by the 
operator 'pW'pi'^)^ jg relevant, while the remaining kinetic or non-hydrodynamic component 
is negligible. Some justifications for this assumption are provided in the discussion section 
of the paper. Then, the second term on the right hand side of Eq. ffTTj) is neglected when 
substituting it into Eq. (!56|) and the solution of this equation can be written down by simple 
inspection, 

0h(ci, C2) ~ Pl'^pf x(ci)<5(ci - C2). (57) 

Use of Eqs. (146|) and (l50l) into Eq. (pT]) gives the fluctuating Navier-Stokes equation for the 
velocity field. 



5u{k, s) H ; — k • 5uj{k, s)k 

a 



W{k,s), (58) 



ds 2 

where k = k/k and the noise term W{k, s) = —ik ■ R(fc, s) has the properties 

{W{k,s))H = 0, (59) 
11 



{Wi{k, s)W,{k\ s'))h = '^5u,-k'G{\s - s'|)A;2 + . (60) 

The quantity G{s) is evaluated approximately in Appendix O The used approximation, 
which is an exact property for the Maxwell model for inelastic gases reads 



A+(c)A,,(c)x(c) ~ A4A,,(c)x(c). (61) 

In the above relation, A"*" is the adjoint operator of A and the eigenvalues A4 is determined 
self-consistently. The result is 

1 + 02 (a) 



with 



G(s) = —^e'", (62) 

= C« + T^V (83) 

1 + a2{a) 

The explicit form of the functions 02 (a) and I (a) is given in Eqs. (lASP and (ICOp . respectively. 

It is interesting to evaluate the shear viscosity rj given in Eq. flSTl) in the same approxi- 
mation as employed to calculate the function G{s), i.e. by using Eq. (16T]) . In the long time 
limit s ^ 1 it is found in Appendix O that 

1 + a2{a) 



V - , \ - Co • (64) 



This result is very close to the one reported in refs. 2l| and 22|, both being practically 



indistinguishable for a ^ 0.6. Moreover, Eq. f lMl) coincides with the expression derived 



m 



25| by using a modified Sonine expansion in which the gaussian is replaced by the 



(approximated) distribution of the HCS, x(c)- 



V. LANGEVIN EQUATION FOR THE ENERGY FIELD 

First consider the term 6(^{k,s) given in Eq. fl27|) . It will be evaluated here by making 
an approximation in the same spirit as Eq. (I6T]) . namely by writing 

A+{c)Uc)^X3Uc) = -jUc). (65) 
Then, by realizing that Eq. ( 1271) is equivalent to 

6ak,s) = -2{^,\A6F), (66) 
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6C{k,s)^^[5eik,s) + 6pik,s)]. (67) 



it follows that 

The same result is obtained if, instead of Eq. (1651) . the approximation 

{^,\A6F) ^ (^,\AV6F) (68) 

is employed. Therefore contributions to 6({k,s) from V±6F are being neglected in Eq. 
( 167|) . This includes, in particular, terms proportional to the gradients of the fluctuating 
hydrodynamic fields (see, for instance, Eq. f j4T]) ). As for Eq. ( !6T|) . the approximation given 
)y^Eq. f l65p becomes an exact relation for the inelastic Maxwell model for granular gases 



24j |. Moreover, it has been shown that the linear transport coefficients associated with the 



gradient expansion of the average cooling rate are very small as compared with the similar 



contributions coming from the hydrodynamic fluxes [2l| . Actually, the only contribution to 
the hydrodynamic equations from the average cooling rate that is kept in practically all the 
literature is the one of zeroth order in the gradients. 

Next, the heat flux term 6(p defined in Eq. f l24l) has to be evaluated. It contains the 
function S(c) that verifies a relationship similar to Eq. (jHj), 

J dc^{c)^is{c)=0, (69) 

for all the hydrodynamic modes ^^{c). Therefore, Eq. f l24|) is equivalent to 

6(p{k,s) = J dc'E{c)V±SF{k,c,s), (70) 

and by means of Eq. ( H3l) . the energy flux can be decomposed as 

S(l){k,s) = Si(f){k,s) + Z{k,s), (71) 

where 

6i(j){k, s) = ds' J dci:{c)e''^^''\-ik-c)V5F{k,c,s- s'), (72) 

and 

Z(fc, s) = ds' j dcT,{c)U{k, c, s')V±S{k, c, s-s'). (73) 
The above expressions are valid up to first order in k. A direct calculation gives 



V6F{k,c, s) = 6p{k,s) 



6(c) + ^ 



6iv{k,s)-Uc) + ^6e{k,s)Uc), (74) 
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and taking into account the isotropy of the operator A(c), Eq. (172|) can be rewritten as 
5i0(fc,s) = ds' J dc'E{c)e''^^''\-ik- c) 



X |5p(fc, s - s')Uc) + ^ Mk, s-s') + 5p{k, s - s')] ^3(c)| 



(75) 



The next task is to evaluate the functions 6p{k, s — s') and 6e{k, s — s') as functions of the 
fluctuating hydrodynamic fields at time s to lowest (zeroth) order in k, in order to have an 
expression for 5icf){k, s) valid to first order in k. Use of the balance hydrodynamic equations 
(I20])-(I22]) leads to 

5p(fc,s-s') ~5p(fc,s), (76) 
5e{k, s-s') + 6p{k, s-s') ^ e^°''/2 [6e{k, s) + 6p{k, s)] 



r rfsie^«^^/2^,(fc,si). (77) 

J s-s' 



-Co(s-s')/2 



When the above expressions are substituted into Eq. f[75|) . two contributions physically 
rather different are identified. One of them is of hydrodynamic character and similar to the 
expression for the fluctuating heat flux for molecular gases derived by Landau and Lifshitz 



Is 



9[, while the other one is an intrinsic noise term following directly from the inelasticity of 
collisions, 

Si(p{k, s) = 6[^^ct){k, s) + 5f V(fc, s). (78) 
After simple manipulations, the flrst contribution can be expressed as 



5i^V(fc,s) = ds' j dc-E{c)e''^^''\-tk ■ c) 



X i^6p{k, s)ei(c) + ^e^°^'/2 [^^(^.^ + ,)] ^^(c^ 
= —Kik6e{k, s) — {p — k) ik6p{k, s). (79) 

The right hand side of this equation has the same form as the (generalized) Fourier law 
for dilute granular gases 2l|, |22|. It involves two transport coefficients: the (thermal) heat 



conductivity k, and the diffusive heat conductivity p. Their expressions are: 

^(s) = ^ / ^cS(c)-$3(c,s), (80) 

Ks) = \j rfcS(c) -$1(0,5), (81) 
14 



with 



$i(c,s)= rds'e^'^(^)ei(c)c + 2$3(c,s), (82) 

^0 

$3(c,s) = \ r ds' e^'^^^^^^^-/'ks{c)c. (83) 
^ Jo 

As usual, k(s) and Jl{s) are expected to reach steady plateau values for large enough s, when 
the hydrodynamic description is accurate. Both transport coefficients have been evaluated 



in the ffist Sonine approximation 2l|, |22| 



The noise term in Eq. (!78|) is given by 

6[^^cf>{k,s) = ^J^ ds' j t^cS(c)e"'^(")zfc-ce-^«("-"')/2 



X 



r dsi e^«^^/25,(fc,si)6(c). (84) 

Js-s' 



Substitution of Eqs. (|67|), ([n]), (JZHj) and ([79]) into Eq. ([22]) yields 



\ds 2 ) 



6e{k, s) + ^6p{k, s) + i '^—^k ■ 5u:{k, s) 



+ ^ A;2 [K6e{k, s) + (Jl - K)5p{k, s)] = £{k, s), (85) 

(Jj 



where 



£{k, s) = S,{k, s)-'^k- Z{k, s)-'^k- 6[^^(l){k, s) (86) 

is identified as the total noise term. The functions Z{k,s) and 6[^^(f){k,s) are defined in 
Eqs. (173|) and (^^, respectively, and the intrinsic inelastic noise term S^{k, s) is given in Eq. 
f l28|) . It is trivial to verify that 

{£{k,s))H = 0, (87) 

while the corresponding correlation function is calculated in Appendices |D] and [E] under well 
defined and controlled approximations, which will be made explicit also below. The result 

reads 

{£{k,s)£{k\s'))H ^ ^Ua)assia)6{s - ^04,-^' + ^ ^ j S^,-k'k' 



{d + 8) ^ ^ , 2rfCo(a)a33(a) [1 + 2a2(a)] \ ^^^|,_,,| 



X <! 1 + — ^ — a2(a) + ^ y— )■ e 

\-^5\ — Co/ ^ 



(88) 



vahd for s,s' ^ 1. In this expression 033(0) and A5 are given by Eqs. (IDSP and (IDlSp . 



respectively. The two main approximations made to derive the above expression are similar 
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to those leading to Eq. (160|) and (l62l) . Firstly, the non-hydrodynamic component of particle 
velocity correlations is neglected once again. Secondly, it is written that 



A+(c)S(c)x(c) ~ A5S(c)x(c), (89) 

and As is consistently obtained by means of Eq. flD12|) . 

When the expressions of thermal heat conductivity k and the diffusive heat conductivity 
Jl, Eqs. ( IHOj) and ( IHTl) respectively, are evaluated using the two above approximations, the 
results read (see Appendix [D]) 

^ {d + 2) [1 + 2a2{a)] 



2(2|A5|-Co) 



(90) 



^^^_(. + 2)[2_+a,(.)]_ 
4|A5| 

The above expressions hold in the limit of large time s. The above values for the thermal 
transport coefficients are indistinguishable of the results found in the first Sonine approx- 



imation 
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22| for all values of a. On the other hand, they are not equivalent to the 



expressions reported in ref. 25|, although both are very close for a ^ 0.65. 

Finally, to close the description provided by the fluctuating hydrodynamic equations 
derived here, the correlation between the fluctuating force appearing in the velocity equation 
and the fluctuating force in the energy equation is needed. It is verified in Appendix [F] that 

{W{k,s)S{k',s'))H = 0, (92) 

as expected because of symmetry considerations. 



VI. FINAL REMARKS 



The objective here has been to derive fluctuating hydrodynamic equations for a dilute 
granular gas by extending methods which are familiar for normal gases. The main result 
obtained is the set of coupled Langevin-like equations (12(1 . (ISSl) . and (15^ . describing the 
time evolution of the fluctuating hydrodynamic fields. These equations contain two kind of 
terms. There are terms which describe the hydrodynamic part of the fluxes and of the cooling 
rate. The former have the same form as the macroscopic hydrodynamic fluxes, involving 
the Navier-Stokes transport coefficients. They have been made explicit in the equations. 
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On the other hand, there appear the non-hydrodynamic components of the fluxes as well 
as some additional terms in the equation for the energy due to the dissipation in collisions. 
These latter terms have zero average and have been combined all together to define the noise 
terms in the Langevin-like equations. The auto-correlation functions of the noise terms in 
the velocity and energy equations are given in Eqs. (I^U|) and f l88p . respectively, while the 
cross correlation function between both noise terms is indicated to vanish in Eq. (!92|) . 

Generalizing fluctuating hydrodynamics to granular fluids entails several important dif- 
ferences from normal molecular fluids. Primary among them are the following. 

1. The homogeneous reference state about which fluctuations are considered is not the 
Boltzmann equilibrium state, but the HCS. Its distribution function is not a simple 
function of the global invariants. Moreover, this reference state is time-dependent, 
although it is possible to consider a stationary representation of it by using an appro- 
priate time scale (see Eq. ([2])). 

2. Although the noise in the fluctuating inelastic Boltzmann equation is white, i.e. delta 
correlated in time, the noise terms in the equations for the velocity and energy fields 
have finite relaxation times. Moreover, the amplitude of their correlation functions 
is not determined by the Navier-Stokes transport coefficients, but involve new coeffi- 

n 

cients. In other words, the fluctuation-dissipation relations of the second kind |27| are 
not verified. 

3. The noise term in the equation for the energy has a zeroth order in the gradients 
contribution. This noise term is intrinsic to the inelasticity of collisions and has no 
analogue in molecular fluids. 

The expressions for the two-time correlation functions of the noise terms have been com- 
puted using some approximations. It has been assumed that the velocity correlations be- 
tween particles can be accurately approximated by their hydrodynamic part, identified as 
their projection onto the subspace of distributions generated by the hydrodynamic eigen- 
functions of the linearized inelastic Boltzmann collision operator. The second approximation 
used consists in dealing with the hydrodynamic fluxes as if they were left eigenfunctions of 
the above mentioned linear operator. The two hypothesis can be justified on the basis of 
the following features: (1) they have been shown to lead to some predictions that are in 
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very good agreement with molecular dynamics simulation results. This includes the fluctu- 



ations of the total energy of the system 



13 



14| and also of the transversal component of 



the velocity fleld 10|, (2) if the approximations are applied to the formal expressions of the 
Navier-Stokes transport coefficients, very accurate expressions are obtained as discussed in 
previous sections of this paper, and (3) the second approximation mentioned above is an 
exact property in the case of the inelastic Maxwell model kinetic equation. 

Some comments on the context and utility of the results in this work seem appropriate. 
The analysis has focussed on the fluctuations of the hydrodynamic flelds in the HCS. In 
many experimental conditions of interest the system is far from a global homogeneous state. 
Nevertheless, the reference state studied here is relevant locally for more complex and realistic 
conditions, as it is the case of the equilibrium state in molecular fluids. For example, the 
transport coefficients such as the viscosity obtained here are the same functions of density 
and temperature as those in the associated nonlinear hydrodynamic equations applicable 
under more general conditions. Thus the context of relevance of the equations derived in the 
present work are expected to transcend the limitations associated with the state considered 
and extend to states for which the nonlinear Navier-Stokes are required to characterize the 
macroscopic hydrodynamic flelds. 
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Appendix A: Dimensionless linear Boltzmann collision operator 

The dimensionless binary collision operator Tq{ci, C2) for inelastic hard spheres or disks 
is deflned by 

To(ci, C2) = j dS- e(ci2 ■ ct)ci2 ■ [a-Wci, C2) - l] , (Al) 

where C12 = ci — C2, dB is the solid angle element for the unit vector ct, B is the Heaviside 
step function, and 6^^(ci,C2) is an operator changing all the functions of ci and C2 to its 
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right by the same functions of the precolhsonal velocities cl and C2, given by 



* ,-1 l + a.^ 

= b„ Cl = Cl (cr • Ci2)cr, 

za 

* ,-1 l + a 

C2 = b^ C2 = C2 + — (cr ■ Ci2)<T. 

2a 

The expression of the hnearized Boltzmann colhsion operator A(c) is 



(A2) 



A(ci) 



rfC2To(Ci,C2)(l + Pl2)x(c2) 



2 dci 



■ Cl. 



(A3) 



The operator P12 interchanges the labels of particles 1 and 2 of the quantities to its right, x(c) 
is the scaled velocity distribution of the HCS defined in Eq. ( |T2l) . and Co is the dimensionless 
cooling rate for the decay of the temperature of the HCS in the time scale s, 

dTnU) 



Co 



ds 



'1 — a^lvr 2 



<oTh{s) 



dci / dc2c\^x{ci)x{c2)- 



(A4) 



(A5) 



2V{^)d 

Approximated expressions for the distribution function of the HCS and for the cooling rate 
have been obtained by expanding the functions in Sonine polynomials and keeping only the 



fli 



lowest orders [5|, l26i . The distribution function has the form 

x{c) = ^^,[l + a,{a)S^\c-)] 



where 



and 



^^22 8 



a2{a] 



16(1 -a)(l- 2^2) 



(A6) 



(A7) 



(AJ 



9 + 24rf + {8d - 41)a + 30a^ - 30a^ 
Equation ( ]A6p is used all along this paper to carry out explicit calculations. The approximate 
expression for the cooling rate is: 



Co 



727r('^-^)/2(l-a2) 
r (rf/2) d 



1 + 



802(0;) 
16 



(A9) 



19 



Appendix B: Derivation of Eq. ()54p 



From Eq. (HHI) and using Eq. (fT6|) . it follows that for small wavevector k it is 



{R{k,s)R{k',s'))H 




dci j rfc2 A(ci)A(c2) 

s' 



,_,'_s^ + s2)n^nT(ci,c2), 

(Bl) 



where it has been used that to lowest order in k the inhomogeneous linear Boltzmann 
operator A{k, c) can be replaced by the homogeneous one, A(c). The projection operators 
V^l'^ and V^'' act on functions of Ci and C2, respectively. Suppose that s' > s. The integration 
over S2 can be easily carried out to get 



For s' > s ^ 1 and assuming that all the non-hydrodynamic components of r(ci,C2) 
correspond to negative eigenvalues of A(ci) and A(c2), the above relation can be simplified 
to 



(R(fc, s)R{k', s'))h = / dc^ / dc2 A(ci)A(c2)e(^'-^)^('=^)0^,(ci, C2), (B3) 



where 0j|/(ci,C2) obeys Eq. (1561) . The symmetry of the tensor A(c), the isotropy of the 
operator A(c) and the invariance of r(ci,C2) under rotations of Ci and C2 imply that the 
right hand side of Eq. ( 1B3I) must have the form given in Eq. ( |54l) . The case s > s' ^ 1 
follows trivially. 

Appendix C: Approximated evaluation of the function G{s) defined in Eq. (155p 

Using the approximation given in Eq. (1571) into Eq. ( |55ll gives 




(B2) 





(CI) 
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and taking into account the symmetry of the operator A(c) and of the function this 
expression is seen to be equivalent to 

G{s) = j rfcA,.,(c)e^^(^)A,,(c)x(c) 

= (e^^^A,.,x|A,.,x), (C2) 
where A"'"(c) is the adjoint of A(c) defined by 

{g\^h) = {A+g\h)% (C3) 

for arbitrary functions g{c) and h{c) of the Hilbert space. Now the approximation is made 
that Axy{c)x{c) is an eigenfunction of A+, being A4 the eigenvalue, as expressed by Eq. (I6T!) . 
This approximation is prompted by the fact that it is exact for the inelastic Maxwell model 



of granular gases [2J] . Use of Eq. fl6T]) into Eq. ( \C2\i yields 



Gis) = e^^^ I dcclclxic) = e^^^ (C4) 

where the expression of x(c) in the first Sonine approximation, Eq. ( ]A6p has been employed 
to evaluate the velocity integral. 

To determine A4, Eq. (1611) is multiphed by A^y(c) and integrated over c to get 

X^ = —-^—- [ dcc^CyA^{c)A.^y{c)x{c). (C5) 



1 + a2{a) 

The evaluation of the velocity integral on the right hand side of the above expression using 
once again Eq. ( 1A6P is a lengthly but quite standard calculation. The result is given in Eq. 
fl63|l . where 

(2rf + 3-3a)(l + a)7r'¥ 



Ha) 



" ' 23a2(a) 
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(C6) 



2V2d{d + 2)T{d/2) 

In the same approximation as introduced above, the expression for the shear viscosity, rj{s 
given by Eq. (ISTj) in the limit of large s becomes 

^ ^ J^^ds' J rfcA,.,(c)e^'[^(^)-^°/2l6,.(c)c, 
ds' e-'^«/2(e^'^^ (xA.,) |6,.c,) 
(xA.,|e2,x(c)c,) r ds' e-''i^°'^-~^') 



-(xA.,|e2,.c,)(A4-|] . (C7) 
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Taking into account that 

1 



{X^,y\i2,.Cy) = -, (C8) 



Eq. flMl) is obtained. 



Appendix D: Correlation of the noise term in the equation for the energy field 

In the following calculations, it will be assumed without loss of generality that s > s'. 
From Eq. ( I86i) and keeping only contributions up to order k"^ it is obtained that 

{Sik,s)£ik',s'))H ^ {S,{k,s)S,{k',s'))H - ^k- {Z{k,s)Z{k',s'))H ■ k' 



~k ■ {6['U{k, s)S,{k', s'))h - ^k ■ {Z{k, s)S,{k', s'))h. 
a a 



(Dl) 



Upon writing the above equation, it has been taken into account that 6[^^(l){k, s) is at least 
of first order in k, as it is directly realized from its expression in Eq. (184|) . Moreover, it has 
been used that 

{S,ik,s)5(l>{k',s'))H = 0, (D2) 

as a consequence of Eq. ( 1311) . Consider first the self-correlation of the intrinsic noise S^, 
given by Eq. (!30|l or, equivalently, by 

41/2 /■/■__- 
{S,(k,s)S,{k',s'))n = ^Sk,-k'S{s-s') j dcj dc^ Uc)Uc^)V{c,d). (D3) 

The velocity integral appearing on the right hand side of this equation can be evaluated 
exactly using the expression for x(c) in the first Sonine approximation. Nevertheless, for 
the sake of consistency, here the non-hydrodynamic components of the particle velocity 
correlations in the HCS will be neglected, as it was done when solving Eq. ( l56i) . With this 



approximation, it was shown in ref. 



isj that 



{S^ik, s)S,{k', s'))n ^ ^Coa33(a)4,-fc'5(s - s'), (D4) 

with 

,, d+1 d + 2 . . 
033 (a) = + ^^02(0) + 0(a), (D5) 

2 + d - 6rf2 - (10 - 15d + 2d^)a - 2(2 + 7d)a^ + 2(10 - d)a^ 
~ Qd{2d + 1) - 2ci(ll - 2d)a + Uda"^ - 12da^ ' ^ ' 
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Consider next the self- correlation of Z{k, s). To the lowest order in k, it is 



k- {Z{k,s)Z{k\s'))H ■ k' ^ J dc J dc! k ■ i:{c)k' ■ T,{c!) J dsi J ds2 

^gSiA(c)+S2A(c')^^^y^^ C, S - Si)S{k, C', S' - S2))h 



V 



2 



Nd ' 



-Sh-h'k^ / dc / rfc'Sfc) • Sfc') 



X 



/ rfs2e^2^("')e("+"2-"')^(")r(c,c'). (D7) 
Jo 



For s > s' ^ 1, and assuming again that the HCS is stable with respect to homogeneous 
perturbations of the velocity, so that all the non-hydrodynamic eigenvalues of A(c) must be 
negative, carrying out the integration over S2 yields 

k-{Z{k,s)Z{k',s'))H-k' = -^5k,-k'k^ J dc J dc'^{c)-^{c') 

xe(^-^')^(^)0H(c,c'), (D8) 



where (pn is defined in Eq. (]56l) . The right hand side of the above relation can be evaluated 
by means of an approximation scheme similar to that used in Sec. IIVI to compute G{s). 
Since the justification of the approximations to be made is the same as discussed in Sec. IIVI 
it will not be repeated here. To begin with, Eq. (1571) is used to write 

k ■ (Z(fc, s)Z(fc', s'))h ■ k' -Y^^5k,-k'k'' j dc^ic) ■ e(^-^')^('=)S(c)x(c). (D9) 
Now it is assumed that (compare with Eq. (!6T!) ) 

A+(c)S(c)x(c)^A5S(c)x(c), (DIO) 
so that Eq. (1D9P is approximated by 



k-{Z{k,s)Z{k',s'))H-k' -Yls^ _^,k'^e^s-s')x, j dci:\c)x{c 



^]5k,-k'k\d + 2)[2 + {d + 8)a2(«)] e(^-^')^^ (Dll) 



8iV 

To determine the value of A5, Eq. ( IDlOj) is multiplied by Cx and afterwards integrated over 



the velocity c, yielding 

= 

{d + 2)a2{a) 



^5 = , . , I c/cc,A+(c)S,(c)x(c). (D12) 
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From here it is obtained 



- 4J(a) Co(q) , 3Co(a) , . 

{d + 2)a2{a) 02(0;) 2 



J [a) = 



32V2dr {d/2) 

X {16(2 + d){l-a) + a2{a) [70 + 47c/ - 3(34 + 5rf)a]} . (D14) 
The third term on the right hand side of Eq. (IDip is evaluated as follows. It is 
k- {6[^^(f){k,s)S,{k',s'))H = ^J^ dsi J c;cfc-S(c)e('-"i)^(^)ifc-ce"^«^^/2 

X r ds2e'^'''/^S,ik,S2nik',s'))HUc)- (D15) 
Using the approximation defined in Eq. fIDlOp . it is 



/.cE(c)e.-..^<«.C,(c).e<-.^./.cE(c)C3(c) 

= e(^-^^)^-^^ [1 + 2a2(a)] I. (D16) 
Then, use of this result and Eq. (1D4I) into Eq. (ID 1511 yields 

k.{6['^cl.ik,s)Aik',s')), ^ Kd + 2)[l + 2a2ia)]Co{a)aMV'k^ 

x5k,-k'e'''^'^^e''^' / dsie-'^(^^+^°/2). (D17) 
Jo 

In the limit s > s' ^ 1, the above result reduces to 

' ' N - Cn/2) 

The calculation of the last term on the right hand side of Eq. (IDip is much more involved. 
Nevertheless, in the next appendix some evidence is presented yielding the conclusion that it 
can be safely neglected since it leads to a much smaller contribution than the other terms, at 
least for not too strong inelasticity. Then, putting all the results obtained in this appendix 
together, Eq. f l88|) follows. 

The long time limit of the expression for the thermal heat conductivity k, Eq. (IHOj) . is 
easily evaluated within the approximation scheme developed above. 



K 



1 fT , CoV 



~2dr'^f' ldcUcMc)-c. (D19) 
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By using now that 

I dcUcmc) ■ c = 'Mi:^ [1 + 2a2(a)] , (D20) 
Eq. (1D19I) yields Eq. ( l90l) . Equation (I9T1) is derived in a similar way. 

Appendix E: Analysis of the last term on the right hand side of Eq. (ID1|) 

It is 

2? ~ 4? r r 

- —k - {Z{k,s)S,ik',s'))H = -^k- J dsijdcj rfc'c'2S(c)W(fc,c,si) 



xVj_{S{k,c,s-s,)S{k',c',s'))H 



5k,-k'k- J dc J dc' c'^i:{c)U{k,c,s- s')V±T{c,c') 

- jY^2 ■ y j dc' c'^i:{c)U{k,c,s - s') 

xP^[A(c') + A(c)]x(c)5(c-c'). (El) 

In the last transformation, the contribution to V±T{c, c!) due to the velocity correlations in 
the HCS has been neglected, as done everywhere along this paper. Therefore, 

2? ~ 

- -k ■ {Z{k, s)S,{k', s'))h = K,{s- s') + K^{s- s'), (E2) 

with 

Kiis) = ^5k,-k'k- dc dc'c'^^{c)Uik,c,s) 

xP^A(c')x(c)5(c-c'), (E3) 

K2{s) = j^^k-k'k- I dc'E{c)U{k,c,s) 

xV±A{c)x{cy. (E4) 

Using again the approximation in Eq. (165 p . the definition of Ki{s) in Eq. (IE3p can be easily 
rewritten as _ 

Ki{s) - -'^^^6k,-k'k- J dci:{c)U{k,c,s)V^_Uc). (E5) 

Now, to keep up to order fc^, the expansion 

U{k,c,s) ~ e^^^^('=)+ / dsie''^^^^''^V^{-ik-c) 

Jo 

^g(s-si)PxA{c) (Eg) 
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is employed to get 

since the first term on the right hand side of Eq. flE6p gives a vanishing contribution to 
Ki{s) because of symmetry. In principle, the velocity integral on the right hand side of 
the above equality can be evaluated by using the same kind of approximations considered 
along this paper, i.e. by treating T,j.CxX{c) as an eigenfunction of A~^{c) and determining 
the eigenvalue in a self-consistent way. Nevertheless, the contribution given by Ki{s) seems 
to be negligible as compared with the other terms retained in Eq. fl88|) . The argument is 
as follows. Due to the operator V± itself, the time exponential to its right decays with 
eigenvalues corresponding to the kinetic, non-hydrodynamic part of the spectrum of A. 
Because of symmetry considerations, the components associated to the eigenvalue A4 and 
A5 give vanishing contributions to the integral. Then, a sensible estimation of the velocity 
integral seems to be 

J rfcS,.(c)c,Pxe(^-^i)^(^)e3(c) ^ e(^-^^)^'^ J dcT..,{c)c,V^Uc). (E8) 

The underlying assumption is that A5 is an upper-bound for all the kinetic modes. Note 
that, in particular, this is true for A4 as obtained above. In Figs. [1] and |2] the obtained 
approximated expression for Ki{s) is compared with the term proportional to a2(«) on the 
right hand side of Eq. (155]) for a = 0.6 and a = 0.9. Namely, the two plotted quantities are 

A,e'\ A, ^ ('' + 2)(d + 8)a.(a) ^^^^ 

and 

^2se^^% A2 = 2Co j dcJ:,ic)c,V,_Uc), (ElO) 

as functions of s. The latter has been evaluated using the first Sonine approximation for 
x(c). It follows from the figures that for those values of s for which both quantities are not 
negligible, it is |y42(a)|se'^'*'^ ^ |Ai(a)|e'*'^^. It is true that the two functions cross one another 
for large enough values of s, but when this happens both are very small. For instance, for 
a = 0.6 the curves cross at s ~ 6, 7, and their value is of the order of 10~^. For a = 0.9, the 
value of the functions at the intersection roughly — 10~^^. Similar behaviors arre obtained 
for other values of a. The above results indicate that the contribution Ki{s — s') to the 
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0,4 




FIG. 1: (Color online) Comparison of the dimensionless quantities ^le'^^^ (solid blue line) and 
A2se'^^^ (dashed red line) as a function of the dimensionless time s for a = 0.6. The definitions of 
Ai and A2 are given in Eqs. ()E9P and (jElOp and the case d = 2 has been considered. It is seen 
that the contribution Ki{s) associated to the amplitud A2, given in Eq. ()E3p can be accurately 
neglected. 

correlation of the noise term in the energy equation can be safely neglected. The term K2{s) 
given in Eq. (lE4p can be analyzed in a similar way, reaching the same conclusion. This 
justifies neglecting the last term on the right hand side of Eq. (IDip . 

Appendix F: Correlation between the noise terms in the velocity and energy equa- 
tions 

Taking into account that 5[^'*</)(fc, s) is at least of first order in k, it follows from Eqs. 
f HHj) . (!86|) and the definition of W{k,s) given above Eq. (!59l) that to Navier-Stokes order, 
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FIG. 2: (Color online) The same as in Fig. [H but for a = 0.9. 



i.e. to second order in the gradients, it is 



{W{k, s)S{k', s'))h = -ik ■ (R(fc, s)S,{k', s')h - -^k ■ (R(fc, s)Z{k\ s))h ■ k' . 



Consider first 



{Rij{k,s)S,{k',s')H = \ J J dc' c''^Aij{c)U{k,c,si) 



xV^{c){S{k,c,s-si)S{k\c',s'))H 
2V^ 



Nd 



Sk,-k' J dsi j dc J dc c''^Aij{c)U{k, c, Si) 



x6{s - si - s')r(c, c') 



(Fl) 



(F2) 



The right hand side on the above equation contains the integral 



Iij= / dc / dc' c'^Aij{c)l({k,c,si)r{c,c') 



(F3) 
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Since there is already an explicit k factor in Eq. (IFip . this quantity is needed to first order 
in k. Therefore, the approximation 



- ds2 e'^'^^'^^^^'^^V^ik ■ ce^^^-^^)^^^^'^)^^ (F4) 
Jo 



'0 

is used. The zeroth order in k contribution to lij is 



.c/..'.^A«(c)e..-.r(c.c', 
^/.cA,,,.)e..-./.c'c»f(..c',^0. (F5) 

since Tic,d) is invariant under rotations of c and d. The first order in k contribution to 



lij is 
It is 

and also 







l\^ = -ik- I ds2 / rfcA,j(c)e^^^(")cP^(c)e(^i-^^)^('^) j dc' c'^T{c, c'). (F6) 



(^2(c)|e^^(^) I dc'c''T{c,c'))=0 (F7) 



J (icA,,(c)e^^('=)ce^(c) = 0, (F8) 

for P = 1,3. Therefore, the operator V±{c) on the right hand side of Eq. flF6l) can be 
omitted, and the expression is seen to vanish since it has the form 

= -tk- j^' ds2 j dc A,,(c)e^2^(^)c(7(|c|, Si - S2) = 0. (F9) 

To compute the second term on the right hand side of Eq. (]F1|) to Navier- Stokes (fc^) 
order, the correlation function appearing there must be evaluated to order zero. Then the 
projection operators can be eliminated by using the same kind or arguments as above and 
it is obtained 

y2 rs r-s' 

{Rij{k,s)Zi{k',s))H ^ -TT^k,^k' / dsi / ds2 6{s - si - s' + S2) 

Jo Jo 

^ J dcj rfc'A,j(c)Sz(c')e^i^^'=)e^^^('^')f(c,c') = 0, (FIO) 



29 



again because of symmetry considerations as a consequence of the isotropy of A and the 
invariance of r(c, c') under rotations. This completes that proof that the correlation in Eq. 
(IFip vanishes. 
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